Level Statistics and Localization for Two Interacting Particles in a Random Potential 
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We consider two particles with a local interaction U in a random potential at a scale L\ (the one 
particle localization length). A simplified description is provided by a Gaussian matrix ensemble 
with a preferential basis. We define the symmetry breaking parameter fi oc U~ 2 associated to the 
statistical invariance under change of basis. We show that the Wigner-Dyson rigidity of the energy 
levels is maintained up to an energy _E M . We find that _E M oc when F (the inverse lifetime 

of the states of the preferential basis) is smaller than A2 (the level spacing), and oc 1//J, when 
r > A2. This implies that the two-particle localization length L2 first increases as \U\ before 
eventually behaving as U 2 . 

PACS numbers: 72.15, 73.20 



> 

o 
o 

<N 
O 

vo 

OS 
-I— > 



C 

O 

o 



X 



For a single particle diffusing in a disordered system 
of size L smaller than the one particle localization length 
L±, there are two characteristic energies: the Thouless 
energy E c — hD/L 2 and the level spacing Ai w B\jL d 
(Bi , D and d are the band width, the diffusion constant 
and the system dimension, respectively). If one writes 
the distribution of energy levels as a Gibbs factor of a 
fictituous Coulomb gas, the corresponding pairwise inter- 
action for levels with separation e < E c coincides |lj with 
the logarithmic repulsion characteristic of the matrix en- 
sembles which are statistically invariant under change of 
basis, e. g. the Gaussian Orthogonal Ensemble (GOE). 
For e > E c , the level repulsion vanishes more or less 
quickly, depending on the system dimension. The dimen- 
sionless conductance g\ is given by E c /Ai, This ratio is 
the single relevant parameter in the scaling theory of lo- 
calization. In quasi-one dimension, the size where g\ ss 1 
defines L\. In three dimensions, the mobility edge is 
characterized by g\ ss g c where g c is of order 1. 

We shall generalize those concepts to two particles with 
a local (repulsive or attractive) interaction. This two in- 
teracting particle (TIP) problem has received a partic- 
ular attention since Shepelyansky B pointed out that 
certain TIP states may extend over a scale L2 much 
larger than L\. Shepelyansky's original reasoning con- 
sists in mapping the problem for L L\ onto a random 
band matrix model with a superimposed diagonal matrix 
(SBRM- model). Imry || used later the Thouless scaling 
block picture to arrive at precisely the same results as 
Shepelyansky. The smearing due to the interaction of 
the energy levels within L\ was estimated using Fermi's 
golden rule, yielding L2 oc U 2 . This derealization effect 
has been confirmed by transfer matrix studies ^,^], and 
unambiguously illustrated from numerical studies || of 
rings threaded by an ^4i?-flux. However, in one dimen- 
sion, for system sizes which can be numerically inves- 
tigated, one obtains M L 2 oc \U\ contrary to Fermi's 
golden rule, and a disorder dependence Q Li oc L" with 
a w 1.5 — 1.7 and not 2, as predicted by Shepelyansky 



and Imry. 

To understand those contradictory results, we study 
the TIP energy level statistics at a scale L\ in order to 
identify the energy which plays the role of E c in this case, 
and to determine its dependence on U . For the original 
TIP-problem, we assume a tight-binding model (2) on a 
d-dimensional lattice (Lf sites p where the random poten- 
tial is taken with a box distribution of width 2W). The 
nearest neighbor hopping term takes a constant value 
V = 1 and U is the on-site interaction. Assuming two 
electrons with opposite spins, we consider only the sym- 
metric states. The TIP-Hamiltonian || can be written in 
the basis of the N = Lf(Lf + l)/2 (symmetrized) prod- 
ucts of one particle states \AB). We denote by R p a the 
value on site p of the one particle eigenstate with en- 
ergy ca- In this basis, the diagonal terms are dominated 
by one particle contributions + eb and the interaction 
Hamiltonian yields a full matrix (for L < L{) with entries 
U ■ Qaba'B' = U J2 P R p aR p bR p a' R p b' ■ The magnitude 
of those terms is of order U/L^ d ^ 2 with a random sign 
which does not preserve the sign of the interaction. 

Before considering the TIP-Hamiltonian, it is instruc- 
tive to discuss a simplified matrix model where the corre- 
lations between matrix elements are neglected: an ensem- 
ble of real symmetric matrices G with independent en- 
tries, characterized by Gaussian distributions with vari- 
ances < G% Sf/3 (B 1 = Wd + 2W) and < G 2 V] >« 
U 2 /L\ d for the diagonal and off-diagonal terms, respec- 
tively. The averages are set to zero, which neglects a shift 
of the diagonal terms by an amount U ■ Qabab ~ U /Lf 
assumed to be much smaller than B\. These shifts pre- 
serve the sign of U and for large U, eventually split the 
energy band into two parts. For the sake of simplicity, 
we ignore them, restricting us to small U and to a Gaus- 
sian matrix with preferential basis (GMPB)-model which 
has been used previously 0] to study the GOE to Poisson 
crossover for the level statistics, and to define a maximum 
entropy model j^] where the range of the level interaction 
depends on a parameter. When < G 2 i > is very large 
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FIG. 1. The strength function for a TIP-Hamiltonian (5 
by 5 lattice in the metallic regime W = 2, V = 1). Di- 
amonds, squares and triangles are for U= 0.2, 0.6 and 1.0, 
respectively. Lines are Breit-Wigner distribution functions, 
fitted to the numerical data. The inset shows how F depends 
on U. The line represents F = U 2 /22. 

compared to < G?- > , one has indeed a strongly prefer- 
ential basis and it is convenient to re-order the diagonal 
terms such that G\\ < G 22 < ■■■ < Gnn- Those Gu 
may be considered as the positions of the energy lev- 
els in the zeroth approximation, when the interaction 
with the other states is neglected. The small coupling 
terms GV, spread those basis states over T/A 2 neighbors. 
A 2 ps \Jl-K < G\ >/N is the level spacing and T can be 
estimated using Fermi's golden rule: r pa 2tt < G?- > 
/A 2 . Diagonalizing G by an orthogonal transformation 
O (G = OGdO 1 , where Gd is a diagonal matrix with 
real entries E a ), we consider the strength function intro- 
duced by Wigner Pw (E,p) = £^ =1 2 p J{E - E a ). If 
one averages over the ensemble, one gets for the eigen- 
vector amplitudes < 2 a >= A^ 1 < pw(E a ,p) >. For 
L ^> L\, pw(E,n) has been found ]9|-[Tl|] in agreement 
with the Breit-Wigner form 



< pw{E, n) >- 



2tt[(E - G n 



I*/4] 



(1) 



We have checked that this holds too for L < L\, and we 
show in Fig. 1 that this Breit-Wigner form of the eigen- 
states characterizes also the original TIP-Hamiltonian, 
once the shifts U ■ Q nn of the quasi-energies H nn are 
taken into account. This shows us that a basis state 
\n) = \AB) (i. e. an eigenstate of the TIP-Hamiltonian 
for U — 0) becomes delocalized over r/A 2 of its neigh- 
bors (i.e. over the basis states \n') where H n > n > is close 
to H nn = eA + £b), with a Lorentzian shape centered in 
H nn + U ■ Qnn- T plays the role of a localization length 
in the preferential basis, and is given by Fermi's golden 
rule for small values of U (U < 2 in Fig. 1). 

Having understood how the eigenstates are delocalized 
by the interaction over the preferential basis (see also 
Refs. HQ), we focus our attention on the energy levels. 



We introduce a symmetry breaking parameter fi in the 
probability density 



p(G) oc exp 




(2) 



with cr 2 ps Bf/3 and cr 2 /(2(l + //)) ps U 2 /L\ d . When 
p = 0, one recovers the GOE-ensemble with pgoe(G) oc 
cxp(— tr(G 2 )/2a 2 ). When p ^ 0, there is a factor p M (G) 
which removes the statistical invariance under change 
of basis. Expressed Q in eigenvalue-eigenvector coor- 
dinates, it reads 
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The question is to understand how this additional fac- 
tor, after integration over the matrices O (distributed 
with Haar measure p(dO) over the orthogonal group) can 
destroy the logarithmic level repulsion coming from the 
measure p{dG) = n« </3 \E a - E \ n„ dE a p(dO). This 
will allow us to identify the characteristic scale E„ below 
which one recovers the GOE rigidity, and above which 
the levels become uncorrelated. Two cases have to be 
considered: 

(i) r < A 2 . The Gy are so small that one can just con- 
sider the coupling between two nearest neigbor diagonal 
entries: i. e. a 2 x 2 matrix which can be diagonalized 
by a rotation of an angle 0. One finds || J ddp fl (G) = 
f(x) = exp(— x) ■ Io(x) where x — p ■ e 2 /(8er 2 ), e denot- 
ing the separation of the two coupled levels. For x < 1, 
f(x) ps 1 and decreases as l/\/x for x > 1. This gives 



E n _ ^/8a 2 /p 



x 



N 



(4) 



For e < E n , one has the GOE statistics, while for e > E^, 
the levels are uncorrelated. 

(ii) r > A 2 . Many neighboring Gu are coupled by the 
off-diagonal terms. First, we consider the case where 
e = \E a — Ep\ < T: i. e. the case where the two 
corresponding eigenvectors have a strong overlap. As- 
suming that the eigenvectors \O a ) have non-zero co- 
ordinates of order 2 na ps A 2 /r over T/A 2 neighbor- 
ing basis states only, one gets 2 va 2 v p ~ A 2 /T, 
and exp (—pe 2 A 2 /(2cr 2 r)) ps 1, independent on e(< T). 
Writing O — exp A, with A a real antisymmetric ma- 
trix (p(dO) — n a< pdA a/3 ), one can see that the small 
fluctuations of the A na around their typical values will 
not yield a correction to the GOE level repulsion. This 
means that there is no coupling between eigenvalues and 
eigenvectors as far as e < T 



E n with 



E» N 2 
A 2 p 



(5) 
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FIG. 2. E 2 for the GMPB-ensemble (N = 500). Di- 
amonds, squares and triangles are for fi= 30000, 5000 and 
1000, respectively. The inset shows how depends on [i. 
The triangles give the energy where E 2 is 20 percent above 
the GOE value. The solid and the dotted line represent 
E fi /A 2 = 0.039iV 2 /^ and £ , M /A 2 = 0.19N/^/Jl, respectively. 
The squares result from a fit E 2 = (E/E^) a valid for 
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FIG. 3. E 2 for the symmetric states of a 2d TIP-Hamil- 
tonian (ring with 10x10 sites, W = V = 1). The ring is 
threaded by a magnetic flux $ = <&o/4. Diamonds, squares 
and triangles are for U= 0.25, 0.75 and 2.0, respectively. The 
inset shows how E\j depends on U. The data are obtained as 
described in the caption of Fig. ^| The dotted line (solid line) 
corresponds to E v /A 2 = 1MU/V (E v /A 2 = 2.1(U/V)' 2 ). 



now, instead of N/y/Jl previously. When |e| > T, the 
eigenvectors do not overlap and the levels should become 
uncorrelated. In Ref. ||, it was noted that if O pa w 

5 p , a +A pa where A pa < 1, J2p=i °la°l/3 ~ 2A lf3> which 
gives a l/|e| factor, after integration over A a p. This level 
attraction exactly compensates the level repulsion due to 
fi(dG). Qualitatively, one can adapt this reasoning to 
produce the requested level attraction, after integration 
over the eigenvectors. Quantitatively, the calculation of 
the exact form of the level interaction as a function of e, 
taking into account the Breit-Wigner form of the eigen- 
vectors, is postponed to a further study. 

We have carried out a numerical study of the GMPB- 
ensemble (N = 500) as a function of fi, to illustrate 
the two regimes. The number variance £2 (£7) (vari- 
ance of the number of levels in an energy interval E) is 
shown in Fig. 2. For small energy intervals, H^iE) coin- 
cides with the GOE-logarithmic increase observed when 
fi = 0. For larger energy intervals, ^(E) can be fitted 
by (E/E^W, which gives a first method for calculating 
Ep. A second method consists in calculating the energy 
interval where £2 is above the GOE-curve by a certain 
threshold (e. g. 20 percent). Note that those methods 
give a non zero E^ (depending on the chosen threshold) 
even for uncorrelated levels, which has been subtracted 
from the data. The inset of Fig. 2 confirms that the 
two methods are in agreement and exhibits the predicted 
crossover for £/, when r « A2 (/i « -/V 2 ), from a iV 2 //i- 
dependence (small fx) towards a N/ y^-dependence. The 
/i-dependence of the exponent a (see inset in Fig. 4) de- 
pends on the exact form of the pairwise level repulsion. 

We now turn our attention to the TIP-Hamiltonian in 
two dimensions. The corresponding £2 is shown in Fig. 3, 
for different U. We have obtained the same curves 



for U = 1 and U = — 1, and consider in more detail 
repulsive interactions. We took rings containing 10 x 10 
sites threaded by a magnetic flux, so that the level statis- 
tics should have a GUE behavior for energy intervals 
E < E v = E^, with (1 + ll)- 1 w W 2 l{L\ d B\). Ex- 
cept this trivial change from orthogonal to unitary sym- 
metry, the similarity with the curves calculated for the 
GMPB-model is very striking, if one disregards large en- 
ergy intervals and the [/-dependence of the exponent a 
(see inset of Fig. 4). The first difference can be explained 
from the fact that in the TIP Hamiltonian, the one parti- 
cle level rigidity cannot be ignored when E > Ai (one has 
A1/A2 superimposed GUE series when U = 0), correla- 
tions which are neglected in the GMPB-model. Another 
difference results from the shifts U ■ Q n n of the diagonal 
terms which become important when U is large. The 
crossover value U c w V2(8Vd + 4W)/(^ttL^ 2 ) between 
the two regimes is of order 1, for the considered parame- 
ters. When U > U c , we observe the J7 2 -behavior of E\j. 

A similar study in one dimension is very instructive. 
When W = V = 1, we have L\ w 25, which gives again 
U c ~ 1. As expected, one can see in Fig. 4 that E\j cx \U\ 
when \U\ < U c , but when U > U c , the splitting of the 
energy band occurs, and E\j decreases. Note that one 
recovers a Poisson statistics when U is very large (for 
d = 2, there is only a saturation of E\j). For d = 1, this 
means that one can couple only two basis states within 
L\ w 25, with a small enough value of U to justify the 
simplified GMPB ensemble. The observation of the U 2 
behavior of E\j requires larger values of L\ in d = 1 than 
considered in the numerical studies H-lfl . 
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FIG. 4. Exs for a Id TIP-Hamiltonian (W = V = 1 
and L = Li = 25 sites). The data are obtained as de- 
scribed in the caption of Q. The dotted line corresponds to 
Exj/IS.2 = 1.2U/V. The inset gives the exponent a occurring 
in E2 oc E a at large E. Squares are for the Id TIP Hamilto- 
nian, plotted as a function of UVL/V. Diamonds represent a 
for the GMPB ensemble (see Fig. 0), as a function of N/y/Ji. 

We now follow the argument developed by limy || to 
estimate the localization length L 2 . First, we consider a 
series of building blocks of size Lf. Ttj is the smearing of 
the TIP levels of one of the blocks, due to the interaction- 
induced coupling with the neighboring block. For such a 
quasi-ld wire, the dimensionless conductance at scale L\ 
is given by 

The factor 1/2 gives || the right limit when ?7 — *■ and 
A is a constant. Obviously, one should have Ttj = E\j. 
When T\j > A 2 , T\j is given by Fermi's golden rule, the 
case considered in Ref. H and we only discuss the case 
Tu < A 2 , where Tu ~ \/U 2 / L\ d . Physically, this means 
that U is so small that it couples only a single TIP state in 
one of the blocks to another TIP state in the next block, 
giving rise to "Rabi oscillations" between those two cou- 
pled states. The inverse life time is no longer given by the 
square of the coupling term, as in Fermi's golden rule, but 
by its absolute value. In addition, we have shown that 
this inverse life time gives the scale below which one has 
a GOE spectral rigidity, extending the known results for 
non-interacting particles to the TIP problem. 

We conclude by discussing a few implications for TIP 
localization. For d = 1 and U < U c (U c ~ 1 when 
W « V « 1, see Fig. 4), one gets {L 2 /L x ) w 1/2 + 
A(\U\I B\)\fL\ which is in agreement with the depen- 
dence on U observed in the numerical studies ||. The 
conjecture proposed in Ref. || gives L\ instead of \[L\. 
As noted in Ref. |Q, the distribution of the QabA'B' is 
far from being Gaussian, which can matter as far as the 
description of the Li-dependence by the GMPB-model is 



concerned. However L 2 oc L x ' is close to the behavior 
observed in Ref ||. The [/-dependence is not affected 
by this consideration and is correctly described by the 
GMPB-model. 

If one considers two quasi-particles above a Fermi sea, 
one should replace § in Eq. (§) A 2 by A 2 (E) w Af/E 
where E is the total excitation energy. One immediately 
obtains that the quasi-particle conductance g^E^Li) is 
of order of g 2 (Li) when E sa B\, which gives L\[E ss 
B\) = L2, in agreement with Ref. [p"3f . Similarly, in three 
dimensions, Imry's relation (E m2 ~ (Bi/\U\)E^/ 2 ) be- 
tween the one quasi-particle mobility edge E m i and the 
two quasi-particle mobility edge E m2 does not change 
when U < U c [v denotes the critical exponent associated 
with Li). 

In summary, we have shown that the basic concepts de- 
veloped for non interacting particles can be naturally ex- 
tended to M = 2 interacting particles, after the changes 
Ec — > E^ and Ai — * A 2 . A similar conclusion has been 
obtained from a non linear a model description of the 
TIP-Hamiltonian Jlij, when L > L\. Moreover, our ap- 
proach can be easily extended to an arbitrary number M 
of particles. 
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